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ABSTRACT 


We measure the two-point clustering of spectroscopically confirmed quasars from the 
final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales 
of 4 < s < 22 /i -1 Mpc. The sample covers 6950 deg 2 (~ 19 (/i _1 Gpc) 3 ) and, over the 
redshift range 2.2 ^ z ^ 2.8, contains 55,826 homogeneously selected quasars, which is 
twice as many as in any similar work. We deduce 5q = 3.54 ± 0.10; the most precise mea¬ 
surement of quasar bias to date at these redshifts. This corresponds to a host halo mass of 
~ 2 x 10 12 h~ 1 M@ with an implied quasar duty cycle of ~ 1 percent. The real-space 
projected correlation function is well-fit by a power law of index 2 and correlation length 
ro = (8.12 ± 0.22) /i _1 Mpc over scales of 4 < r p < 25 /i -1 Mpc. To better study the evo¬ 
lution of quasar clustering at moderate redshift, we extend the redshift range of our study to 
z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 ^ z ^ 3.4. 
We find no significant evolution of tq or bias over this range, implying that the host halo mass 
of quasars decreases somewhat with increasing redshift. We find quasar clustering remains 
similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is 
monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with pre¬ 
vious BOSS measurements, but they yield more precise constraints based upon a larger and 
more uniform data set. 


1 INTRODUCTION 

It has long been known that a supermassive central bl ack h ole oc¬ 
cupies every massive galaxy (e.g.. lKormendv & Richstonell 19951 ). 
and that the gal axy and the black hole properties are strongly 
correlated (e.g., iMagorrian et alj 1 19981 : iFerrarese & Merritt] l200d : 

iGebhardt et all l2000i ). Galaxies and their black holes therefore 

appear to have co-evolved. Although, there are other non-causal 
expla nations that have been invoked to explain these relations 
(e.g., iKormendv & Hoi 120131) . This situation could arise be¬ 
cause feedback from quasars partially regulates star formation, 
and thus the properties of galaxy bulges, through merger-driven 
winds, gas accretion or stochastic mergers (e.g., iNandra et alj 


1 20071 : Isilverman et alj 120081 : IShankail l2009h . Alternatively, this 
relationship could occur because both supermassive black holes 
and their host galaxies are governed by the prop erties of 
their parent dark matter halo (e.g., iFerraresel 120021) . Quasar 
populations show a high level of clustering, and it is in¬ 
ferred that they occupy dark matter ha los of ~ 10 12 h~ 1 Mq ) 


at most epochs ( e. g.. iPor ciani. Magliocchetti & Norberg 


da Angela et al. 12008 :IRqss et alj2009l:1white et a 1 1 2012 ). Star for 


Croom^eAlJJ^^^^^rs^et^,. 2006, jTOOTckJjihenjt al 
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mation is also most efficient in halos of this mass (e.g., Moste r et alj 
20ld:lBehroozi, Conroy & Wechslerll2010l : lBethermin et alj|2012l : 
Viero et alj|2013bl) . providing further circumstantial evidence for 
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a link. Further study of the relation between quasar activity and 
galaxy evolution is clearly warranted. 

The geometry of the cosmological world model is now 
known with remarkable precision, and even systematic dif¬ 
feren c es are at the few percent lev el (e.g,, HinshawetalJ 
20131: IPlanck Collaboration et al.1 120141 : lAnderson et al. 1 20141 : 


Aubourg et alJ 20141 ). The spectrum of primordial fluctuations and 


its evolution over cosmic history is also becoming increasingly 
well constrained. Armed with this knowledge, it has become in¬ 
creasingly meaningful to study galaxy formation over a range 
of redshifts—particularly at high redshift, where subtle changes 
in the cosmological model can greatly affect predictions about 
galaxy evolution. Driven by this better cosmological understand¬ 
ing, coupled with large surveys of active galaxies over a range of 
wavelengths, supermassive black holes are increasingly being fit 
into a wider cosmological framework for galaxy formation (see 
lAlexander & Flic kox 2012, for a comprehensive review). 

From the quasar perspective, a wider understanding of 
the interplay between galaxies, star formation and black 
holes has been studied extensively in extraga lactic surveys 
(again see Alexander&HictoJ 120121 as well as Hogkmsetd] 


20071: IConrov & Whitel 120131: IVeale. White & Conrovl I20l4 


Caplar, Lilly & Trakhtenbrotl |2014|) . In addition to direct mea¬ 


surements of the relationship between the luminosity of Active 
Galactic Nuclei (he nceforth AGN) lumin o sity and star formation 
in g a laxies (e.g.. lAlexander & Hickoxl 1201 & iMullanev et alJ 
12013: iRoss et al.ll2012bl) .~ great progress has been driven by two 
statistical measures of the quasar population: the luminosity 
function and clustering. The quasar luminosity function is increas¬ 
ingly well-understood (e.g., see IRoss et alJ 12013d . for a recent 
study) and there are now a large number of quasar clustering 
measurement s in th e o ptical at 2 < 2.2 acros s a ran ge of scales 


(e.g.. IPorc iani. Ma gliocchetti & Norb er3l2004 : Croom etakll 20051: 


IPorciani & Nor herd 20061; iHennawi et alJ I 2 OO 6 I: iMvers et all 


l2007al lbL 120081 : Ross et al] 1 20091) . Together, these measurements 
provide an increasingly coherent characterization of the dark 
matter environments inhabited by quasars and what fraction of 
quasars ignite in tho s e environments (i.e,, t h e quasar duty cycle; 
Cole & Raised 198^: Haiman & Huil 2001: Martini & Weinberg! 


200 ll: I White. Martini & Cohnl 120081: IShankar et alJ l2010h . Mea- 


surements of quasar clustering in the optical have informed 
similar studies at 2 < 2 


24.3 h x Mpc at 2 > 3.5. It remains unclear exactly how these 
large correlation l engths transition to rp ~ 5/t~ 1 Mpc (7 ~ 2.0) 
at 2 < 2 (e.g.. IPorciani & Norbergl 1200a ). In addition, quasar 
clustering measurements in samples selected beyond the optical 
are now being made at moderate redshift, and comparisons be¬ 
tween these samples and (rest-frame) UV-luminous quasars have 
historically been critical to interpret the quasar phenomenon at 
2 < 2. Star formation has been shown to occur most efficiently 
in halos of similar mass to those occupied by quasars at 2 < 
2 (~ 10 12 fc^Mo), and as measurements of star formation as 
a function o f dark matter halo mass are now being pushed to 
2 > 2 (e.g., IViero et ail 1 201 .3 al . 1201 4l) it will be intriguing to see 
whether this trend continues to higher redshift. Finally, models 
of how quasars populate dark matter halos, for instance, empir¬ 
ical approach es such as the halo occupation distribution frame- 
work (HOD; Seljak[2000|^ Peacock^&Sjmitl^ 2000|^ Bensonetal] 


200(j^ SroccimaiToet^y|2001^ White 1 Hemauist & SnringefeOOll : 

iBerlind & Weinberdl 2002 l : ICoorav & Shethl I2OO2L and references 

therein) have been successful i n modeling galaxies a t low 
redshift (e.g ., Zheng et aD 120051 : IZheng. Coil & Zehavil 120071 : 
IZehavi et all 1201 lh . but are now being updated to model 


Degraf et al. 201 la; Chatterjee et al. 

2013; Richardson et al 

White et al. 2012; Conroy & White 

2013; Richardson et al. 

Veale, White & Conrov 2014; Caplar, Lilly & Trakhtenbrot 


2013 


It is not clear to what extent such HOD models describe rare 
and heavily biased quasars at high redshift, and it is critical to 
push observations to 2 ~ 3 as the most sophisticated hydrody¬ 
namic simulations being used to interpret the HOD are now being 
pushed from high red shift down to 2 < 4 (e.g., DeGr afet all 20121 : 
iDi Matteo et all2012l) . 

Given the importance of new observational constraints on 
quasar clustering at moderate redshift, we have embarked on mea¬ 
surements of quasar clustering at 2 < 2 < 3 using a set of 
q uasars fr om the Baryon Oscillation Spectroscopic Survey (BOSS; 
iDawson et alj2013j) part of the third phase of the Sloan Digital Sky 
Survey (SDSS-III: lEisenstein et al.ll201 ll) . Although the main goal 
of BOSS was to constrain cosmological models, such a large sample 
of quasars offers the chance to study galaxy evolution through mea¬ 
surements of quasar clustering in the largely unsampled redshift 
range of 2.2 < 2 < 2.8. BOSS reaches more than tw o magn itudes 


(e.g., Allevato et al. 

201 ~i|; Krumpe et alJ 2012 j)_in the radio 

2002; Schneider et al. 2013). This increased depth, combined with 

(e.g., Shen et al. 

2009 

). in the infrared (e.g., Donoso et al. 20141; 

improved targeting algorithms, allowed BOSS to spectroscopically 

DiPompeo et al. 

2014allbl; Hickox et al. 2011) and, in general. 

confirm 15 times as many 2 ~ 2.5 quasars as earlier iterations of 


in AGN _samples selected across the electromagnetic spectrum 
(e.g.. iHickox et al]l2009l) . Beyond 2 < 2 relatively weak direct 
constraints on quasar clustering at 2 > 2.9 also exist dShen et alJ 
120071 120 I 0 I ) and similar high-redshift measurements are starting to 
be made for quasar samples identified outside of the optical (e.g., 
lAllevato et all2014h . 

It is now critical to understand quasar clustering in the range 
2 < 2 < 3 (or at “moderate redshift” for quasars) for a number 
of reasons. Luminous quasa rs peak in number density at moder¬ 
ate redshift (e.g., Ro ss et al] 1 201 3al ). so this is the most promis¬ 
ing epoch for understanding how quasar feedback and dark mat¬ 
ter environment are related. Measurements of quasar clustering are 
generally sparse at moderate redshift, and such redshifts may be 
key to linking the relatively weak constraints at 2 > 3 to the 
more accurate measurements at 2 < 2. Most notably, IShen et ail 
J2007I) found (for a power-law index of 7 = 2 . 0 ) that quasars 
cluster with r 0 = 16.9/i _1 Mpc at 2 ~ 3 growing to ro = 


the SDSS. 

I White et alJ d 2012 l) reported preliminary quasar clustering 
measurements from BOSS using the autocorrelation of ~ 27,000 
quasars spread over ~ 3600 deg 2 . IWhite et al. |H3) found a 
quasar bias of 5 q = 3.8 ± 0.3 at 2 ~ 2.4. Using a similar 
BOSS sample, Font-Ribera et al. (2013) cross-correlated ~ 60,000 
quasars with the Ly a forest and derived &q = 3.64+g|g at 
2 ~ 2.4. Font-Ribera et al. could use a larger sample size be¬ 
cause their adopted cross-correlation technique is insensitive to the 
angular select ion function of BOSS quasars. In this paper, we up¬ 
date the Iwhite et~aD I 2 OI 2 I ) BOSS quasar autocorrelation results 
using the considerably larger and final BOSS sample. BOSS has 
now spectroscopically confirmed ~ 160,000 2 > 2.2 quasars 
to a depth of g < 22 (and/or r < 21.85). About half of these 
quasars are uniformly-selected, with about two-thirds of that uni¬ 
form, or “CORE” sample populating a large contiguous area of 
~ 7500 deg 2 in the North Galactic Cap (NGC: IRoss et alj[2012bl ). 















































































































































































Clustering of intermediate redshift quasars 3 


This ~ 7500 deg 2 is by far the largest contiguous area ever used 
for 2 > 2.2 quasar clustering measurements, and represents a con¬ 
tiguous area ~ 3x larger than used in IWhite et alJ d2012h . The 
main BOSS quasar sample and different subsamples that we use 
for our measurements are described in lj2] The clustering method 
and uncertainty estimation are presented in fj3]and the clustering 
results are presented in S|4] The evolution of real and redshift-space 
correlation functions in terms of redshift and luminosity are dis¬ 
cussed in Sj5] In Sj6]we interpret our measurements in terms of char¬ 
acteristic halo masses and associated duty cycles, with particular 
attention to the observed (lack of) luminosity dependent cluster¬ 
ing. We summarize our conclusions in Sj7] Appendix A introduces 
the software that we used to generate a uniform random catalog 
from the survey mask. We adopt a ACDM cosmological model 
with fl m = 0.274, Qa = 0.726, h = 0.7 and as = 0.8 as as- 
sumed in other BOSS analyses (e.g., Whiteet aT]|20 1 ll : iReid et all 
l2012l:IWhite et all2012l:)Anderson et all2012h . We express magni- 
tudes in the AB system l Oke & Gunnll 1 983lh 


2 DATA 

2.1 BOSS Quasars 


The SDSS-l/ll/lII imaging surveys mapped over 14,000 deg 2 of 


w 

sky using the SDSS_ camera jGunn et al. 1998!) mounted on a 2.5- 
meter telescope dGunn et al.l2006l) located at the Apache Point Ob- 
servat ory. Of the unique imaging area released as part of SDSS 
DR8 dAihara et al.]|201lh . the BOSS survey targeted quasars over 
~ 7400 deg 2 of a region around the North Galactic Cap (NGC) and 
~ 2600 deg 2 of a region around the South Galactic Ca p (SGC ). The 
B ary o n Oscillation Spectroscopic Survey (BOSS: iDawson et alJ 
120131) . one of four SDSS-III surveys dEisenstein et al.l20Tll) . has the 
primary goal of measuring the clustering of luminous red galaxies 
and neutral hydrogen to determine the cosmic distance scaltf]. As 
part of this effort, BOSS measured spectroscopic redshifts of 1.5 
million galaxies and 190,000 quasars using twin multi-object fiber 
spectrographs dSmee et al.1120131) . In this paper, we focus on using 
these data to measure quasar clustering near z ~ 2.5. 

BOSS quasars are selected from point sources in 5-band ugriz 
imaging (Fukug ita et al .1 1 99fil) that meet a selection of photometric 
flag cuts (see the appendixes of lBovv et alj201 lhlRoss et alj2012bl) 
and that fall within magnitude limits of i > 17.8 and dereddened 
g SJ 22.0 or r ^ 21.85. Targeting quasars near z ~ 2.5 is diffi¬ 
cult, as metal-poor stars typically have colours very close to those 
of unobscured quasa rs around redshift 2.7 (e.g.. iFan et al.lfl999l : 
iRichards et alJ Eoolh . and quasars at redshifts of 2.2 to 2.6 can 
have similar colours to quasars at redshift of about 0.5 because host 
galaxy light mimics t he Lyman-o Forest entering the u-band (e.g .. 


S: 


iBudavari et alJl200ll : IRichards et al.ll200ll : I Weinstein et al.ll200- 
These colour similarities are exacerbated near the faint limits of 
imaging, where flux uncertainties are larger. To circumvent these 


Bow et all J201 ll) developed the XDQSO algorithm to tar¬ 


get quasars for BOSS spectroscopic observation. BOSS targeted all 
sources with a (2.2 < z < 3.5) XDQSO probability of pqsoraidz 
> 0.424 tuned to produce 20deg -2 quasar targets over the BOSS 
footprint dRoss et alj|2012bl) . This was known as the “CORE” sam¬ 
ple. 

The primary mission of the BOSS quasar survey was to study 


quasar clustering in the Lyman-a Forest via methods that are in¬ 
sensitive to the quasar selection technique. As such, a heteroge¬ 
neous range of multi-wavelength data was used to target BOSS 
quasars, and a number of methods were incorporated into the sur¬ 
vey in addition to XDQ SO dRichardset alj |2009l : lYeche et al .11201(1 : 
iKirkpatrick et alJ[201 ll : see lRoss et alJl2012bl) . This heterogeneous 
or “BONUS” approach augmented the CORE sample. 

In this paper we focus on BOSS quasars spectroscopically 
confirmed prior to 5th April 2014, which is essentially the final 
BOSS quasar compilation. These quasars correspond to BOSS spec¬ 
tral reduction pipeline version 5-7-0 (spall-v5_7_0). Figure [7] 
shows the angular distribution of our quasar sample. The yellow 
background shown in Figure |T| represents the survey area, and 
the point sources are the CORE quasars over the redshift range 

2.2 < z < 2.8. 


2.2 Targeting Completeness 

Quasar clustering measurements are sensitive to the angular com¬ 
pleteness of the quasar sample (the “mask”). Quasar spectra in 
BOSS were obtained using a set of spectroscopic tiles that can be 
mode l ed as polygons using the MANGLE softwa re dBlanton et alJ 
120031 : iTegmark et alJl2004l : ISwanson et al Jl2008ll . The fraction of 
targets (quasars) in an area of the sky that has been covered by a 
unique set of these tiles (called a s ector) is used to define the com¬ 
pleteness of the sample. Following lWhite et al.ld2012l) . we limit the 
survey to sectors with targeting completeness ( f comp ) greater than 
75%. About 60% of such sectors are 100% complete, and the mean 
area-weighted completeness of the sectors we use is ~ 95%. A veto 
mask is applied to remove the regions in which quasars can not be 
observed such as around bright stars, and near the centerposts of 
the spectroscopic plates dWhite et ali2012l) . 

Note that fiber collisions make it impossible to obtain spec- 
tra for quasar pairs closer than 62" ex cept in overlapping regions 
dBlanton et al.l2003l:lwhite et alj2012h . At the mean redshift of our 
survey (z ~ 2.43), 62" corresponds to a transverse comoving dis¬ 
tance of 1.2 (i _1 Mpc. We have not corrected for this fiber collision 
effect, as it is significant only near the lower edge of the range of 
scales we consider; this is a region where the errors from shot noise 
are already quite large. 

We adopt the spectroscopic redshift reported by the spec¬ 
tral reduction pipeline spAll-v5_7_0 to compute the comoving 
distances. To avoid using any target with a problematic redshift, 
we only accept quasars with zWARNING=0, indicating a confident 
spectroscopic classification and redshift measurement for quasar 
targets dBolton et al.ll2012l f 2 l 

Unless otherwise specified in the text, we also limit our clus¬ 
tering measurements to just the NGC region of Figure[T| as the NGC 
represents the largest contiguous area in BOSS. Excluding the SGC 
quasars from the original overall CORE sample of 69,977 quasars 
significantly affected the clustering signal both in real and redshift 
space. Previous analyses have also found unexplained differences 
between clustering measurements in the BOSS NGC and SGC re¬ 
gions (e.g. JWhite et alJl2012D . The target selection papers for the 
coming eBOSS survey (Myers et al. 2015, Prakash et al. 2015, in 
preparation) suggest that the SGC contains areas in which true target 
density is difficult to regress, which may explain such discrepan¬ 
cies. TableQ]documents how our available sample sizes are affected 


https://www.sdss3.org/surveys/boss.php 


1 


http://www.sdss3.Org/dr9/spectro/caveats.php#zstatus 
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Figure 1. Aitoff projection of the angular distribution of quasars in the final BOSS sample in equatorial coordinates. The background (yellow) depicts all sectors 
of the BOSS survey, as well as a random catalog (for all survey sectors) produced using the BOSSQSOMASK software (see the Appendix).The (blue) dots are 
BOSS CORE quasars in surve y sector s that have a completeness greater than 75%. The CORE sample corresponds to objects that have been uniformly-selected 
using the XDQSO algorithm feovv et al.l201 lh . 



Figure 2. Normalized redshift distributions for the 55,826 NGC-CORE, 
30,551 NGC-BONUS and 86,377 NGC-CORE+BONUS quasars over the red- 
shift range 2.2 ^ z ^ 2.8. Also indicated is the redshift distribution gener¬ 
ated for the random catalog that is used to mimic NGC-CORE quasars. 


by restricting to the NGC, to completenesses of f co mp > 0.75, and 
to ZWARNING=0. 

2.3 Division by Redshift and Luminosity 

To investigate the evolution of quasar clustering, we extended the 
redshift range by relaxing the upper limit from 2.8 to 3.4 and split 
the NGC-CORE sample into three redshift bins that are populated 
by nearly the same number of quasars (see Table[4]l. We also check 
whether the clustering signal is luminosity dependent by dividing 
the NGC-CORE sample into three luminosity subsamples based on 
the absolu te magnitude of the quasars in the i-band. Following 
iRoss et aD d201 3al> . we k-c orrect all of our quasar magnitudes to 
z — 2. Figure[3]presents the absolute magnitude distribution of the 



2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 


Figure 3. The absolute magnitude distribution for BOSS quasars. The upper 
panel displays the 10 th , 25 th , 50 th , 75 th and 90 th percentiles of Mj vs. 
redshift. The right-hand side histogram for the upper panel is the normal¬ 
ized distribution of Mi for the full NGC-CORE sample over 2.2 ^ 2 : ^ 3.4. 
The middle panel shows the same percentiles for Mi — M* ? i vs. redshift, 
where is the characteristic luminosity derived from Eqn.[7] Similar 
to the upper panel, the right-hand side histogram for the middle panel is 
the normalized distribution of Mi — for the full NGC-CORE sample 
over 2.2 ^ 2 : ^ 3.4. The three vertical dashed lines depict the lower limit, 
median and the upper limit for the full redshfit range of the samples stud¬ 
ied in this paper (2.2, 2.49 and 3.4 respectively). The lower panel is the 
normalized redshift distribution of quasars over the range 2.2 ^ 2 : ^ 3.4. 
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Properties 

# of quasars 

XDQSO targets: 

In the CORE sample 

107455 

In the BONUS sample 

104971 

Targets with zWARNING=0: 

In the CORE sample 

103552 

In the BONUS sample 

89044 

Targets with f co mp > 0.75 and zWARNING=0: 

In the CORE sample 

97644 

In the BONUS sample 

79699 

NGC^CORE sample 

78151 

NGC BONUS sample 

60788 


/comp > 0-75, 2.2 2.8: 


In the CORE sample 

70892 

In the BONUS sample 

44786 

/com P > 0.75, 2.2 ^ z ^ 2.8, zWARNING=0: 


In the CORE sample 

69977 

In the BONUS sample 

41019 

NGC CORE sample 

55826 

NGC BONUS sample 

30551 

/comp > 0.75, 2.2 ^ z ^ 3.4: 


In the CORE sample 

94306 

In the BONUS sample 

60277 

/com P > 0.75, 2.2 < z sg 3.4, zWARNING=0: 


In the CORE sample 

92347 

In the BONUS sample 

54614 

NGC CORE sample 

73884 

NGC BONUS sample 

40781 


Table 1. Number of objects in the SDSS DR12 BOSS XDQSO sample with 
different cuts on redshift and completeness. 


quasars in the final BOSS NGC-CORE sample with respect to the 
characteristic luminosity of quasars at every redshift in this range: 

= -21.61 - 2.5 (kiz + k 2 z 2 ) - 0.71 , (1) 

where k\ = 1.39 and k 2 = —0.29 for z ^ 3 and fci = 1.22 
and k 2 = —0.23 for z ^ 3 (from Icroom et all |2004 modi¬ 
fied by ICroton et al.l [2006 ). Quasars at higher redshift are fainter, 
and fainter quasars have larger flux uncertainties which translate 
into imprecision in measured colours. So, at higher redshift the 
^correction becomes increasingly uncertain, (e.g. iRichards et alj 
120061 ) . 


3 CLUSTERING MEASUREMENTS AND ERROR 
ANALYSIS 

To determine how much the distribution of quasars deviates from a 
homogeneous distribution, it is compared to a “random”, or control, 
set of objects. This random catalog has all of the characteristics of 
the data, including the selection function, except for its clustering. 
Figure [2] shows how the redshift distribution of the random cata¬ 
log we have produced mimics that of the final CORE BOSS quasar 
sample. Appendix A introduces the code that we use to make the 
random catalog from the survey mask and quasar redshift distribu¬ 
tion. 


The two-point correlation function £(*) is defined as the prob¬ 
ability, in excess of random, of finding a pair of objects at a sepa¬ 
ration x : SP = nS VJ 1 + fix'll JP eeblesIl 1 98(1) . We use the esti¬ 
mator introduced by lLandv & Szalav 1 1993 ) to calculate the real 
and redshift-space two-point correlation functions (2PCFs), £(a;), 
of quasars in our sample: 


^ (DD(x)) - 2 (DR(x)) + (RR(x)) 

i[X) ~ (RR(x)) ’ { ) 

where DD(x) is the number of data-data pairs with separation be¬ 
tween x and x+Ax and DR(x) is the number of data-random pairs 
also separated by x. The random catalog is constructed to be larger 
than the data catalog to reduce Poisson noise due to the pair counts 
that include random points. Consequently, there are normalization 
factors of the form -/Vr/IVd that scale the counts appropriately if 
the random set has a different size than the data set. The angled 
brackets denote the suitably averaged pair counts. Throughout this 
paper we use random catalogs 12 to 20 times larger than our quasar 
samples. Tests with random catalogs up to 30 times larger than the 
data show little statistical improvement in the clustering signal, but 
increase the computing time for our analyses. 

Using the comoving distances to each of the two objects in a 
pair (si and S2), and the angular separation between them ( 6 ), we 
define each pair's comoving separation along (7r = |si — S 2 I) and 
across (r p = [si + s 2 \9/2) the line of sighjj. We thus derive the 
two-dimensional correlation function £(r p , 7r). Since the effects of 
redshift space distortions are limited to the line of sight (radial) 
direction, we can integrate £(r p , n) over tt in order to eliminate the 
redshift-space distortion effect: 


cut 

w p( r p) = 2 / £(7r,r P )d7r, (3) 

Jo 

where the upper limit, 7r cu t, is the distance at which the effects 
of redshift-space distortions become negligibl^f]- We refer to w p as 
the projected correlation function. Averaging over all pairs whose 
total redshift-space separation lies in a bin produces the monopole 
moment of the redshift-space correlation function, £(s). 

The size, bias and number density of the sample together 
with the redshift range over which the clustering signal is being 
measured, can change the optimum value of 7r cu t. For instance, 
70fe~ 1 Mpc was used for the 2dF QSO survey of 23,3 38 quasars 
1 daAngela et alj l2005t) . ~ 100/i -1 Mpc was used by IShen et alj 
1 20071) for 6,109 quasars at z 5s 2.9, and 50/t _1 Mpc was used 
for the sample of 27.129 DR10 quasars in the SDSS-III/BOSS sur¬ 
vey i jWhite et alj l2012h . The value of 7r cu t should be chosen to 
balance the advantage of integrating out redshift-space distortions 
against the disadvantage of introducing noise from uncorrelated 
line-of-sight structure. Moreover, integrating over a wider range of 
7r means projecting a larger 3D space into smaller r p bins (e.g., 
IWhite et al]2012l) . After testing a wide enough range of upper lim¬ 
its to truncate the integral, we found 7r cu t = 50 /i -1 Mpc to be the 
optimum limit beyond which the £ (r p , 7r) values become negligible 
for our sample. 

We estimate the un certainties using inver s e-variance-weighted 
jackknife^ resampling dScranton et alj 120021 ; IZehavi et alj 1 20021 : 
iMvers etaDl2005h . We divide our sample into N pixels on the 


4 Also denoted by a or R in the literature. 

5 The quantity w p (r p ) integr ates out the effects of both peculiar veloc- 
ities and redshift errors (e.g., ICroom et all 120051 ; Ida Angela et alj [20051 : 
I White et all 20 12b . 
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sky, then create N subsamples by neglecting pixels one-by-one 
and using the remaining area to calculate the 2PCF. We pixelize 
our sample into angular regions that are specified by HEALPIX0 
pixels dGorski et al 112005 ) with N s ^ e = 4 (~ 15° on a side). To 
compensate for missing areas, we merge pixels with fewer random 
points than two-thirds of the meatQ The inverse-variance-weighted 
covariance matrix elements, Cij , are 


Cij=C( n, rj )=J2 

L 


RR,L(ri) 


^ RR(n) 

RR L {rj) 


RR(r. 


[cj L (ri) - w(ri)] 

[UL^j) - ( 4 ) 


where the L subscript refers to the neglected pixel and the random 
pair ratio, RRl / RR, compensates for the relative area through the 
relative random catalog size in each pixel (e.g lMvers et al.ll2007al ). 
One-standard-deviation jackknife errors, oy, are the square root of 
the diagonal elements of the covariance matrix. We assume the er¬ 
rors are Gaussian, and thus compute the goodness-of-fit or likeli¬ 
hood of any model that fitted to the data through 


X = ^2 [w(n) - Wmodel(fi)] Cij 1 [w(rj) - tt> mo dal(»j)] • 
id 

(5) 


4 CLUSTERING RESULTS 

4.1 Real-space correlation function 

To measure the projected auto-correlation function w p (r p ) we use 
a sample of 55,826 NGC-CORE quasars within the redshift range 

2.2 ^ z ^ 2.8 (see TableQ}. 

For sufficiently larg e 7r max , w p (r p ) can be related to the real- 
space clustering as JDavis & Peebles!19831) 



As a first model, we use a simple power law of the form (r/ro) 1 
for the real-space correlation function. For this power-law form, 
Eqn.[6]reduces to 

w P (r P ) _ V¥T[(7- l)/2] fro_y 

r P r[ 7 /2] U; ' 

We approach the fitting process in two separate ways: 

(1) First, we perform a two-parameter fit using the x 2 mini¬ 
mization method to determine both the slope and the intercept of 
the fitted line. We use the resulting slope (7 ~ 1.85) as a fixed pa¬ 
rameter and perform a one-parameter minimization, obtaining a 
scale length of ro = 7.7 ± 0.2 h _1 Mpc with Xied = 0-36 over 9 
degrees of freedom. 

(2) Second, we fix 7 = 2 and fit a 1-parameter model. This 
allows us to evaluate how much clustering amplitude is sensitive to 
the chosen slope of the model. This fit results in a scale length of 
ro = 8.12 ± 0.22 /i^ 1 Mpc with x? e a = 0-25 over 10 degrees of 
freedom. 

Since a slight change in the slope does not make a significant 
difference to the result of the fit for the clustering amplitude, we 


6 http://healpix.sourceforge.net/ 

' Although this process could, in principle, be performed multiple times, 
we only merged the pixels once. 



Figure 4. The projected auto-correlation function for the CORE sample over 

the redshift range 2.2 . z ■< 2.8. The dashed line is the best-fit one- 

parameter power law over the range 4 < r p < 25 /i -1 Mpc. Because the 
two-parameter power-law fit does not have a significantly different slope to 
the one-parameter power-law over the chosen fitting rang e, we adopt 7 = 2 
throughout this paper for consistency with Whi te et a l. <[2012). 


conclude that the clustering amplitude for our quasar sample in the 
fitted range is not strongly sensitive to the slope of the real-space 
correlation function. As we find that 7 = 2 is an acceptable value 
of the slope, we proceed by using 7 = 2 throughout the rest of this 
paper so that our fitted cor relation lengths can be compared more 
easily to other works (e.g lShen et alj|2007l : lwhite et al.ll2012l) . 

Fig-Spresents the results of the fits to the measured w p . The 
fitting range in both approaches is 4 < r p < 25 /i -1 Mpc. We find 
that x 2 maintains an acceptable level, and that our fitting results 
are consistent within the derived uncertainties, for all fitting ranges 
from 4 < r p < 16fi -1 Mpc to 2 < r p < 50/i _1 Mpc. An im¬ 
proved mask and/or deeper imaging for the initial target selection 
of rare and faint high-redshift quasars would be needed to push the 
fitting scales further in order to measure subtle large-scale features, 
but the quasar bias measurement should be well characterized over 
the moderate scales of < 25 fi -1 Mpc. 

Using the method described in tj3] we estimate the covari¬ 
ance matrix of the auto-correlation function (ACF) for NGC-CORE 
quasars. The w p result for the ACF of the NGC-CORE quasar sam¬ 
ple over the fitting range 4 < r p < 25 fi -1 Mpc and its uncertainty 
(a) are reported in the first three rows of Table|2] The remainder of 
the table lists the correlation coefficients as estimated from the co- 
variance matrix calculated via inverse-variance-weighted jackknife 
resampling. 


4.2 Using BOSS BONUS quasars 

As mentioned in tj 2 ] BOSS also surveyed a heterogeneously selected 
quasar sample with a complicated angular selection function—the 
BONUS sample. Although modeling the mask for the BONUS sam¬ 
ple is beyond the scope of this paper, we can cross-correlate a 
BONUS+CORE sample against the CORE sample, using the ran¬ 
dom catalog constructed for the CORE sample, in an attempt to 
improve our measurement precision. Expanding the main sample 
by adding quasars from the BONUS sample might be expected to 
improve the clustering signal (as there are more total pair counts 
in a CORE+BONUS-to-CORE cross correlation than in a CORE-to- 
CORE ACF) but it is worth determining whether such an approach 
is justified by improved statistics. 
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r p 

4.36 

5.18 

6.15 

7.31 

8.69 

10.33 

12.28 

14.59 

17.34 

20.61 

24.50 

Wp 

54.64 

48.87 

43.41 

38.09 

29.77 

26.47 

22.74 

16.68 

11.90 

15.27 

7.80 

a 

12.45 

9.45 

9.56 

8.64 

5.68 

5.88 

4.17 

3.44 

3.16 

3.01 

3.09 


4.36 

1.000 

-0.256 

-0.052 

-0.343 

0.027 

0.485 

0.221 

0.171 

0.174 

0.553 

0.543 

5.18 

- 

1.000 

0.023 

0.152 

0.455 

-0.215 

0.320 

-0.298 

-0.299 

-0.245 

-0.241 

6.15 

- 

- 

1.000 

0.232 

-0.119 

-0.379 

-0.083 

-0.225 

-0.261 

-0.419 

-0.273 

7.31 

- 

- 

- 

1.000 

0.083 

-0.410 

0.095 

-0.050 

-0.958 

-0.606 

-0.579 

8.69 

- 

- 

- 

- 

1.000 

0.472 

0.202 

0.153 

0.656 

0.538 

0.529 

10.33 

- 

- 

- 

- 

- 

1.000 

0.484 

0.382 

0.388 

0.138 

0.119 

12.28 

- 

- 

- 

- 

- 

- 

1.000 

0.031 

0.035 

-0.166 

0.154 

14.59 

- 

- 

- 

- 

- 

- 

- 

1.000 

0.004 

0.323 

0.309 

17.34 

- 

- 

- 

- 

- 

- 

- 

- 

1.000 

0.320 

0.306 

20.61 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1.000 

0.083 

24.50 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1.000 


Table 2. Correlation coefficients as estimated from the covariance matrix computed by jackknife error estimation for NGC-CORE quasars. The first three rows 
are the comoving separation, r p (in /r —1 Mpc), the auto-correlation w p and its uncertainty. The remainder of the table lists the correlation coefficients as 
estimated from the covariance matrix computed by jackknife error estimation. 


Fig- m shows the cross-correlation of CORE+BONUS (86,377 
quasars) against CORE (55,826 quasars) for objects in the NGC 
with zWARNING=0 and completeness greater than 75%. For the 
cross-correlation estimator, the “DD” pairs in Eqn. [2] consist of 
one member from the CORE and another from the CORE+BONUS 
sample. The random catalogue, here, was made for the uniformly 
selected (i.e. CORE) sample, so the “2DR” pair breaks into a “DR” 
pair for the CORE and a “DR” pair from the CORE+BONUS sample, 
where ”R” always denotes the random catalog constructed for the 
CORE sample. In other words, Eqn.[2]for the cross-correlation es¬ 
timator that we adopt is £(*) = (( DcDcb{x )) — ( DcR(x )) — 
DcbR(x) + (RR(x))) / (RR(x)) , where the subscripts “c” and 
“CB” refer to the CORE and CORE+BONUS samples respectively. 
The appropriate data/random normalization factors are included 
here (as in the description of Eqn. [2]). 

The difference between the ACF and cross-correlation at every 
scale is well within the 1 <t uncertainty. However, the cross correla¬ 
tion has slightly better precision on smaller scales, which would 
be a benefit to measurements of. e.g., the one-halo term in the 
HOD formalism. The consistency of the clustering replacing the 
ACF with a cross correlation is a useful check on our analyses, and 
may also represent one way forward for performing clustering mea¬ 
surements using heterogeneous samples. However, given the only 
minor improvement, we continue using the ACF throughout this 
paper. 


4.3 Redshift-space 2-point correlation function 

Using redshift to infer distance for the 2PCF causes a complication; 
peculiar velocitie s introd uce redshift-space distortions in £ along 
the line of sight lKaisei]|l987l ). This effect has utility, though, as 
measuring these redshift space distortions can constrain cosmolog¬ 
ical parameters. 

Fig. [6] shows the redshift-space correlation function £(s) of 
55,826 NGC-CORE quasars over the redshift range 2.2 ^ z ^ 2.8. 
As described in a jackknife resampling is used to derive both 
the uncertainties and the full covariance matrix. The black curve 



Figure 5. Projected cross-correlation between the 55,826 CORE and 86,377 
CORE+BONUS quasars over the redshift range 2.2 $1 z $1 2.8 (red open 
circles). The cross-correlation result is consistent with the auto-correlation 
of the CORE sample (green filled circles) over our scales of interest. As little 
additional precision is introduced by using the full CORE+BONUS sample, 
we restrict our analysis to just the auto-correlation of the CORE sample. 


is the real-space correlation function, f (r) of dark matter com¬ 
puted from HALOFIT Ismith et alJI2003I) modified by the relation¬ 
ship between the redshift-space and real-space correlation func¬ 
tions llKaiseil 1987h 

a s ) = (b 2 + l bf + ^a r) , ( 8 ) 

where / = [D m (»] 0 ' 56 is the gravitational growth factor. The 
model is then fit to the measured correlation function for quasars 
by x 2 minimization. The best-fit bias of quasars relative to the un¬ 
derlying dark matter is 6 q = 3.54 ± 0.1 with x? e d = 1-06 for 
7 degrees of freedom using the full covariance matrix. The lower 
panel of Fig.|6]presents the bias residuals and the horizontal dashed 
line represents the best-fit value of the bias. The vertical dashed 
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r p 

(/i -1 Mpc) 

Wp 

ACF 

a 

Wp 

CCF 

a 

s 

(h~ 1 Mpc) 

£0) 

a 

2.09 

69.353 

19.294 

73.474 

17.498 

3.06 

2.556 

0.609 

2.53 

60.807 

14.491 

64.523 

13.311 

3.70 

1.663 

0.384 

3.06 

57.125 

10.707 

63.607 

10.438 

4.48 

2.406 

0.318 

3.70 

52.908 

8.654 

58.216 

8.133 

5.41 

1.654 

0.213 

4.47 

45.454 

5.981 

46.928 

5.365 

6.55 

1.020 

0.132 

5.41 

39.291 

10.607 

36.604 

11.831 

7.92 

1.034 

0.101 

6.55 

34.392 

9.142 

33.677 

9.015 

9.59 

0.676 

0.069 

7.92 

25.846 

6.051 

28.909 

6.567 

11.60 

0.562 

0.050 

9.58 

23.907 

5.210 

27.066 

5.061 

14.03 

0.433 

0.036 

11.60 

16.056 

3.002 

19.692 

3.129 

16.97 

0.291 

0.026 

14.03 

11.771 

1.854 

16.066 

3.196 

20.54 

0.193 

0.018 

16.97 

13.753 

5.348 

18.433 

8.045 

24.84 

0.150 

0.014 

20.54 

9.404 

3.130 

13.768 

5.087 

30.06 

0.090 

0.010 

24.85 

6.824 

0.958 

9.8349 

3.872 

36.37 

0.069 

0.007 

30.06 

6.554 

2.099 

10.018 

3.256 

43.10 

0.040 

0.006 

36.37 

7.341 

1.422 

10.170 

2.184 

53.23 

0.020 

0.004 

44.00 

4.192 

1.061 

8.122 

5.365 

64.40 

0.010 

0.003 

53.23 

4.483 

1.076 

8.535 

3.208 

77.92 

- 0.000 

0.002 

64.40 

2.576 

0.159 

6.711 

3.004 

94.27 

0.004 

0.002 


Table 3. Projected and redshift-space correlation functions. ACF w p represents the auto-correlation of the NGC-CORE sample, CCF w p is the cross-correlation 
between the NGC-CORE+NGC-BONUS and NGC-CORE samples and £(s) is the redshift-space correlation for NGC-CORE quasars over the redshift range 
2.2 2 . 8 . 



2 4 6 8 10 20 40 60 80 100 

s (h _1 Mpc) 


Figure 6. The redshift-space correlation function £(s) for 55,826 B OSS NGC-CO RE quasars over the redshift range 2.2 ; 2 A) 2.8. The black curve is the 
real-space correlation function, §(r) of dark matter from HAI.OMT dSmith et al.ll2003 I modified by the correction for redshift space distortions from Kaiser 
(1987). The resulting best-fit bias (green curve) is 6 q = 3.54. The lower panel depicts the bias residuals compared to the fiducial fit. The (red) vertical dashed 
lines at ~ 4 h ~ 1 Mpc and ~ 22 h ~ 1 Mpc depict the scale over which we fit the data. 


line indicates the edge of the last bin that is included in the fit. This 
figure demonstrates that scale-independent bias combined with a 
ACDM cosmology is a sufficient fit to our data over scales of at 
least 4 < s < 22 /i -1 Mpc. 

Table [3] presents our measurements of the projected auto¬ 
correlation function and redshift-space auto-correlation function, 
as well as the projected cross correlation with BONUS quasars (see 
S3- The reported uncertainties for all correlation functions are 
the diagonal elements of the covariance matrix derived by inverse- 
variance-weighted jackknife resampling. Our result for the redshift- 
space correlation function is compared to other works in Fig. [7] The 
figure shows some consistency between previous results, although 


our work is in better agreement with a ACDM cosmology on large 
scales. 

In general, the excess of the data above the model seen by 
ourselves and other authors for scales greater than ~ 25 /i -1 Mpc 
is likely to be due to the difficulty in building a random cat¬ 
alog that fully mimics complex and subtle deviations in target 
density that manifest on large scales where the corre lation am- 
plitude is sma l l (e.g . jRoss et alJ 1201 ll : iHo et al.ll2015l) . Notably, 
lAgarwal et alJ d2014l) demonstrate that unknown systematics can 
affect quasar clustering amplitude on very large scales. The dis¬ 
agreement between our adopted model and the data on scales of 
s > 27 /i _1 Mpc is likely to be due to the very small error bars 
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10 100 
s [h _1 Mpc] 


Figure 7. Comparison of our red shift-space correlation function with ear¬ 
lier works . [white et all ll2012h and Karag iannis, Shanks & Ros^ fopfl) both 
used samples of ~ 25,000 CORE quasars from BOS S over the same re dshift 
range we consider in this paper (2. 2 ^ z ^ 2.8). Ishen et all i2007i) used 
~4,400 quasars from SDSS Data Release 5 over 2.9 ^ z < 5.4. The black 
curve is as for Fig. [6] 


on these scales being driven by correlations in imaging systemat- 
ics, which, when jackknifed, lead to an underestimate of the true 
errors. Weighting the target density using systematics maps to see 
if the clustering fit can be pushed to larger scales is currently under 
investigation, but is beyond the scope of this paper. CDM is known 
to be a good description of underlying dark matter and, therefore, 
in this work, we focus on deriving the quasar clustering amplitude 
and slope for the scales on which CDM represents a good fit, and 
on which we therefore believe our data to be uncontaminated by 
subtle systematics. 


4.4 2-D Redshift-space correlation function, £(r p , - k ) 


As mentioned in iJ3] one can project the redshift space separation 
between two objects into orthogonal components along (7rj and 
across (r p ) the line of sight. The elongation of virialized overden¬ 
sities such as clusters of galaxies along the line of sight, caused 
by small peculiar velocities that are not associated with the Hubble 
flow, should be visible in 2-D redshift-space maps as long, nar¬ 
row filaments aligned with the line of sight dTegmark et alj 20041 ). 
To check whether this “Fingers of God” effect is apparent in our 
sample, we measure £(r p ,7r) for our NGC-CORE quasars out to 
scales of ±20fi _1 Mpc, and display the results in Fig. [8] It is 
clear from Fig.[8]that the redshift-space distortions are strongest for 
7r < 20 /i _1 Mpc, confirming that our adopted redshift-space inte¬ 
gration limit of 7r cu t = ±50 /i _1 Mpc is reasonable. Although we 
choose not to further analyze £(r p , 7r) in this paper, Fig.[8]implies 
that studies of redshift-space distortions would be a meaningful av¬ 
enue for further cosmological studies using BOSS quasars. 


20 




1.0 

0.5 

0.0 

-0.5 


- 1.0 


Figure 8. The two-dimensional redshift-space correlation function 
£(r p ,7r) for 55.826 NGC-CORE BOSS quasars. The extension in the line- 
of-sight direction is due to a combination of redshift errors and “Fingers- 
of-God”. 


Az 

z 

# of 

quasars 

# of NGC 

quasars 

n(10 6 ) 
(h -1 Mpc) -3 

2.20 7 z < 2.39 

2.297 

30824 

24667 

4.283 ± 0.027 

2.39 7 z < 2.64 

2.497 

30707 

24493 

3.001 ±0.019 

2.64 7 £ < 3.40 

2.971 

30816 

24724 

1.057 ±0.007 


Table 4. Properties of the three CORE subsamples over the redshift range 

2.2 . z 7 3.4. The columns are redshift range, average redshift and total 

number of quasars for each redshift subsample with or without SOC quasars. 
The final column is the number density of the NGC quasars in each subsam¬ 
ple, together with Poisson errors. 


5 EVOLUTION AND LUMINOSITY DEPENDENCE OF 
QUASAR CLUSTERING 

5.1 Redshift dependence 

In order to have sufficient dynamic range to characterize the evo¬ 
lution of quasar clustering, we extend the redshift upper limit we 
study from 2.8 to 3.4. We divide the resulting sample of 73,884 
NGC-CORE quasar^f] into three subsamples and measure the pro¬ 
jected 2PCF both in real and redshift space for each subsample. 
Table[4]summarizes the properties of the three redshift subsamples. 
The dividing redshifts (2.384 and 2.643) were selected such that 
the three subsamples contain almost the same number of quasars. 

Fig. E shows the result of measuring the real-space 2PCF as 
a function of redshift. Motivated by the fitting result of the NGC- 
CORE sample described in i j4.ll we fit a power law of index 7 = 2 
over scales of 4 < r p < 25 /t _1 Mpc for each of the three redshift 
bins. The scatter of the points (10 points are included in the fit for 
each 2PCF) results in higher uncertainties on the best-fit correlation 
length—but ro remains similar as a function of redshift. 

To evaluate how much fixing the power-law index affects our 


8 All quasars still have zWARNING=0 and completeness greater than 75%. 
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Figure 9. The projected real-space clustering of the NGC-CORE sample in three bins of redshift containing roughly equal numbers of quasars. See Table[4]for 
the sample properties and Table Ofor the fitting results. The red dashed line is the one-parameter fitted power-law with 7 = 2.0 and the blue solid line is the 
one-parameter fitted power-law to the measured correlation function of the main sample plotted here as a reference line. 



Figure 10. The measured redshift-space correlation function £(s) for the three redshift subsamples detailed in Table[4] The fitting range for all three subsamples 
is 3.0 < s < 22.5 h~ 1 Mpc and the fitting results are detailed in Tabled The green curve is the fitted dark matter correlation function to the measured £(s) 
for each sample and the overlaid blue dashed curve is the fitted curve to the £(s) for the main sample, plotted here as a reference for comparison. 


results, we allow 7 to float as a second free parameter in the fit 
over the same range. There is a considerable difference between the 
assumed 7 = 2.0 and the best-fit value for some redshift bins (71 = 
1.57 ± 0.12, 72 = 1.69 ± 0.37 and 73 = 1.88 ± 0.29). However, 
the best-fit values for r 0 remain consistent within the error bars. 

We measure the redshift-space correlation function for the 
same three redshift-based subsamples. Similar to our approach for 
the real-space correlation function, we fit £(s) for the three sub¬ 
samples over the same fitting range as £(s) for the full NGC-CORE 
sample and measure the bias of quasars relative to dark matter in ev¬ 
ery bin of redshift. Fig. [TO] displays £(s) for the three subsamples. 
The model represented by the solid curve is as for Fig. [6] Lower 
number densities for quasars in the subsamples produce noisier co- 
variance matrices. As a result, we only use the diagonal elements of 
those matrices when fitting the real- and redshift-space correlation 
functions. 

The upper part of Table [5] displays the fitting results for the 


evolution of the real and redshift-space correlation functions. In¬ 
cluding full covariance matrices in our fits increases \ 2 without 
producing a meaningful change in the best-fit value for the bias. 
The covariance matrices for £ are less noisy than for w v on larger 
scales—but extending the fits to larger scales for £(s) (beyond the 
vertical dashed lines in Fig. [6] and Fig. 1 10b. still results in a large 
increase in uncertainty for the best-fit bias, presumably because the 
larger scale measurements are not well fit by the same bias factor 
and ACDM model that fits the data at s < 25 ft -1 Mpc. 

As a further test of our overall masking procedures, we cal¬ 
culate the cross-correlation function between quasars in different 
bins of redshift and display the results in Fig. QT] The lack of any 
strong correlation between different redshift subsamples demon¬ 
strates that the pipeline-based redshifts used for this study are of 
sufficient reliability and that the masking procedure we use to pro¬ 
duce the angular selection function for BOSS quasars is robust. 

There have been few measurements of quasar clustering near 
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Figure 11. The one-dimensional projected cross-correlation function w p (x$) between quasars in different redshift subsamples. Lack of any strong correlation 
between redshift subsamples demonstrates that our masking procedure and pipeline-based redshift determinations are robust over our fitting scales. We measure 
only the coordinate transverse to the line-of-sight (r p ), as different redshift bins are correlated along the line-of-sight. The data deviate from the from the dashed 
line (zero correlation) for scales below ~ 1.5 h~ 1 Mpc; we discard these points for all of the fits. This deviation is most likely caused by the fact that BOSS 
can not place fibers at separations smaller than 62" (1.3/i — ^^Mpc at z = 2.5) on a single plate. We note that the slight deviations (< ±0.1) from zero at 
scales of > 1.5 h~ 1 Mpc in the lower two panels are insignificant given the high covariance of the data on different scales. Nevertheless, these deviations 
could indicate low-level, correlated substructure in our masks that might, e.g., depend on how the depth of the quasar target sample changes with redshift. 
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1.0 10.0 100.0 
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Az 

A Mi 

z 

# of quasars 

b Q 

X?ed 

io(h 1 Mpc) 

X r 2 ed 

2.20 ^ 2 ^ 2.80 

-28.74 s: Mi s: -23.78 

2.434 

55826 

3.54 ± 0.10 

1.06 

8.12 ±0.22 

0.25 

2.20 C 2 < 2.384 

-28.70 s:Mj< -23.95 

2.297 

24667 

3.69 ±0.11 

1.55 

8.68 ± 0.35 

0.45 

2.384 ^ 2 < 2.643 

-28.74 sC Mi ^ -24.11 

2.497 

24493 

3.55 ±0.15 

0.61 

8.42 ± 0.54 

0.26 

2.643 C 2 ^ 3.40 

-29.31 ^ Mi ^ -24.40 

2.971 

24724 

3.57 ± 0.09 

0.66 

7.59 ± 0.66 

0.33 

2.20 ^ 2 s; 2.80 

-28.74 ^ Mi < -26.19 

2.456 

18477 

3.69 ±0.10 

2.19 

8.62 ± 0.27 

0.94 

2.20 ^ 2 ^ 2.80 

-26.19 ^ Mi < -25.36 

2.436 

18790 

3.56 ±0.13 

1.71 

7.94 ± 0.41 

1.63 

2.20 ^ 2 ^ 2.80 

-25.36 sS Mi s: -23.78 

2.411 

18559 

3.81 ±0.19 

0.4 

8.29 ± 0.36 

3.29 


Table 5. Clustering results for NGC-CORE sample and subsamples. The first five columns are redshift and absolute magnitude range, their average values 
and the total number of NGC quasars in each of the luminosity and redshift subsamples. Columns 6, 7 and 8 are the bias values that best fit our measured 
redshift-space correlation function, their uncertainty and the reduced x 2 (over 7 degrees of freedom for the main sample and 9 degrees of freedom for the 
redshift and luminosity subsamples). Column 9 is the best fit correlation length which is the result of fitting a one parameter power-law with a fixed index of 
2.0 to the measured real-space correlation function of quasars in each subsample. The last two columns are the uncertainty on the correlation length and its 
reduced x 2 (over 8 degrees of freedom for the redshift and luminosity subsamples and over 9 degrees of freedom for the main sample). The lower limits on 
A Mi for the redshift subsamples and full samples are the highest luminosity quasars that happen to occur in the sample and no imposed upper limit in the 
selection is applied. 


z ~ 2.5, which indicates the difficulties inherent in uniformly- 
selecting large numbers of quasars near this redshift. Fig. ED and 
Fig. El show the reported values of correlation length and quasar 


bias from previous 

works (Porciani. Magliocchetti & Norberg 

2004: Croom et al. 2005: da Angela et al. 

2005i; Mvers et al. 

2006: Porciani & Norberg 2006; Mvers et al 

2007a|; Shen et al. 

2007:; da Angela et alj 

20081; Ross et al. 2009 

White etahl 2012; 

Font-Ribera et al. 2013 

). 


Some caution is required in interpreting Figures El and El 
as different samples probe different luminosity ranges and differ¬ 
ent groups have made different assumptions and choices when in¬ 
ferring ro and bias (e.g., power-law vs. ACDM fit, projected or 


redshift-space fit, range of separations used) and when computing 
the error bars on their results. For Figure 13, we have corrected 
all bias values to our adopted cosmology by multiplying the au¬ 
thors’ reported values and error bars by us,re P {z) / where 

C8,rep(~) is the value of a g at redshift z in the reported cosmol¬ 
ogy and as .Rd(z) is the value for our fiducial cosmology (see §1). 
IShen et alj 1 2007 ) report values of ro for samples in “all fields” and 
“good fields” (with higher photometric data quality but a smaller 
sample) at effective redshifts z = 3.1 (2.9 ^ 2 ^ 3.5) and 2 = 4.0 
(2 ^ 3.5). For corresjtonding bias_values in FigureElwe use those 
calculated bv lWhite. Martini & Cohnl d2008l ) by fitting the Shen et 
al. measurements. 
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A Mi 

Mi 

# of 

quasars 

# of NGC 

quasars 

n (10- 6 ) 
(/i -1 Mpc) -3 

[-28.74,-26.19] 

-26.87 

23809 

18477 

0.982 ± 0.0072 

[—26.19, —25.36] 

-25.77 

23669 

18790 

0.998 ± 0.0073 

[—25.36, —23.78] 

-24.90 

22499 

18559 

0.986 ± 0.0072 


Table 6. Properties of the three NGC-CORE absolute magnitude subsam¬ 
ples over the redshift range 2.2 ^ ^ 2.8. The columns list the absolute 

magnitude range and its average for each luminosity subsample as well as 
the total number of quasars in each subsample with or without SGC quasars 
(3rd and 4th columns respectively). Mi corresponds to the average absolute 
magnitudes of the NGC quasars in each subsample. The last column is the 
number density and the Poisson errors of the quasars in each subsample. 

The BOSS measurements of quasar bias, from this paper and 
from the Lya forest cross-correlation measurement of Font-Ribera 
et al. (2013), are the most precise constraints at any redshift, and by 
far the most precise at 2 > 2. Our value of = 3.54 ± 0.10 for 
the redshift range 2.2 ^ 2 : ^ 2.8 is compatible with Font-Ribera 
et al.’s value of 3.64 ^q 15 at la, a reassuring level of consistency 
given the difference of measurement methods. Our bias and rg val¬ 
ues are higher than most of those measured at 2 < 2, though these 
are typically for lower luminosity thresholds so it is difficult to sep¬ 
arate redshift and luminosity effects. Our redshift subsample mea¬ 
surements are consistent with constant t o or constant bias over the 
redshift range 2.2 < 2 < 3.4. (Because matter clustering grows by 
a factor of 1.35 over this redshift span, constancy of one quantity 
implies evolution of the other, but our measurement errors are too 
large to definitively establish evolution of either.) The SDSS-based 
measurements of Shen et al. (2007) are significantly higher than 
ours. We discuss the marked difference between our measurements 
and the Shen et al. (2007) value at 2 = 3.1 (6q = 6.9 ± 0.7 for 
“all fields”) in Sj6] 


The lack of either redshift or luminosity trends within our sample 
therefore makes the significantly higher clustering of their sample 
(spanning 2.9 ^ 2 ^ 3.5) somewhat surprising. However, we have 
not attempted to measure clustering for a sample that mimics both 
their redshift range and luminosity threshold. 


6 HALO MASSES AND DUTY CYCLES 

Following the approa c h proposed in a range of previous works (e.g. 
Haiman & Huil I 2 OO 1 I : iMartini & Wei liberal 1200 ll : IWvithe & Loebl 
2005h . we now use our clustering measurements to constrain the 
host halo masses of BOSS quasars and the duty cycles of quasar 
activity in these halos. The constraint on halo mass comes from 
the fact that more massiv e halos have higher clustering bias dKaiseil 
11984 Tinker etaklboiol) , so the observed clustering of quasars de¬ 
termines the characteristic mass of their host halos. Moving from 
this halo mass to a duty cycle requires a further assumption to re¬ 
late luminosity and halo mass. Here we assume that this relation 
is monotonic with no scatter—the black hole in a given halo is ei¬ 
ther “on” (i.e. observed as a quasar) or “off’—and we test whether 
this assumption is consistent with our observed (lack of) luminosity 
dependence. 

For each of our samples we compute both a characteristic halo 
mass Mh and a minimum halo mass Mh. min ( see Table [/}. We 
calculate these values using the formula of {Tinker et al.112OI0I) to 
model b(M). Specifically, we use their Eqn. 6 with parameters cho¬ 
sen for an overdensity of A = 200 and with our adopted cosmol¬ 
ogy. We quote a characteristic mass as, sim ply, the mas s that corre¬ 
sponds to our measured bias in the [Tinker et alj l20ld) formalism 
(i.e. bQ = b(Mh )). Alternatively, in our adopted model, Mh,m in 
can be determined from the range of halos that correspond to a 
given bias measurement via 


5.2 Luminosity dependence 

As described in £12.31 we divide the NGC-CORE sample of 55,826 
quasars into three luminosity subsamples (see Table [6j. such that 
there is a distinct difference between the average luminosity for 
each equally sized subsample. Fig. [14] and Fig. |T5] present the re¬ 
sults of measuring the real and redshift-space correlation functions, 
respectively, for these luminosity-based subsamples. As the lumi¬ 
nosity subsamples have even lower number densities than the red¬ 
shift subsamples, the covariance matrices for the luminosity-based 
ACFs are noisier than those for the redshift subsamples, so we fit 
their real- and redshift-space correlation functions using only the 
diagonal elements of their covariance matrices. An assumption of 
diagonal error covariance may be fairly accurate for these sparse 
samples, as shot noise makes a large contribution to the error bud¬ 
get. 

Figure [l6lshows the dependence of ro and b on luminosity for 
these three subsamples, each with an effective redshift 2 ~ 2.4. 
Because of the difference between fitting a power-law (for ro) and 
a ACDM correlation function (for b), the order of points is not the 
same in the two panels, but it is clear that the data are consistent 
with constant ro or b, independent of luminosity over our measured 
range. Solid curves show power-law fits of the form a(L/Lo) rn for 
which we find m = 0.56 ± 0.80 (for ro) or m = —0.027 ± 0.27 
(for b), consistent with m = 0. Dashed horizontal lines show the 
best fit with m forced to zero. Our highest luminosity sample has a 
threshold Mi < —26.19 vs. Mi < —26.6 for Shen et al. (2007). 


b(M > M h , min ) 




(9) 


i.e. our measured bias 6 q = b(M > Mh,m in) is then in¬ 
terpreted to be the mean bias of halos with mass above Mh,m in 
weighted by the halo abundance dn/dM, which we take from the 
[Tinker et ali (2008i ) halo mass function at the mean redshift appro¬ 
priate to each of our samples (see Tabled. 

Associating Mh , m in with the minimum halo mass required to 
host a quasar above the sample luminosity threshold implicitly as¬ 
sumes that all halos above Mh , m in have an equal probability of 
hosting a sample quasar (i.e. constant duty cycle) and that halos 
below Mh,m in can only host quasars below the luminosity thresh¬ 
old (i.e. quasars hosted in halos below a threshold of Mh ,min are 
not present in our sample). 

We determine Mh and Mh ,min and their errors by match¬ 
ing the best-fit values and error range of bQ recorded in Ta¬ 
bleland report these results in Table [7] We have also converted 
the jFont-Riberaet alj d201 3l) measurements at 2 = 2.4 and the 
IShen et all ^2007 4 measurements at 2 = 3.1 and 2 = 4.0 to halo 
masses using the same formalism and display these alongside our 
results in Figure [171 Unsu rprisin gly, our halo mass results are in 
excellent agreement with [Font-Ribera et alj I2013I) . which are also 
drawn from the BOSS survey and which therefore probe a similar 
sample of quasars to that used in our analy sis. I t is more difficult 
to bring our results into accord with IShen et al.1 l2007f) . who mea¬ 
sure host halo masses for 2 ~ 3.1 quasars that are an order of 
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Figure 12. Evolution of the correlation length ro with redshift. We display the results for all NGC-CORE quasars (open orange square) and for the different 
redshift ranges (open black circles) and absolute magnitude bins (filled black triangles) studied in this work (see Table[5j- Results from prior works have been 
recalibrated to reflect our chosen cosmology. The inset provides a more detailed view of the range 2.2 ^ ^ 3.0. 
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Figure 13. Evolution of quasar bias with redshift. We display the results for all NGC-CORE quasars (open orange square) and for the different redshift ranges 
(open black circles) and absolute magnitude bins (filled black triangles) studied in this work (see Table|5j. Results from prior works have been recalibrated to 
reflect our chosen cosmology. The inset provides a more detailed view of the range 2.2 ^ ^ 3.0. 


magnitude larger than our masses at a ~ 3.0 . The observational 
origin of the larger measured masses in lShen et alj l i2007t) is clear 
from Fig.[72]and Fig.|T3] Our clustering bias remains flat with red¬ 
shift and as characteristic halo masses are smaller earlier in cos¬ 
mic history, our flat b(z) implies a dwindling halo mass at higher 
redshift. The biases measured by IShen et all 120071 ). on the other 
hand, greatly exceed our measurements at earlier times, implying 


that b(z) grows steeply with redshift and that quasars would there¬ 
fore occupy higher-mass halos at earlier times. 

The higher biases measured by Ishen et alj 120071 ) ultimately 
have their root in the raw clustering measurements displayed in 
Fig. 0 Much of the difference between our measurement in Fig. 
[7] and that of Ishen et all 120071 ) is driven by data on large scales 
(s > 30/iT 1 Mpc), and, in par ticular, by a single bin around 
~ 35 /t^ 1 Mpc. Ishen et aD 120071) fit their clustering measurement 
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Figure 14. To constrain any luminosity dependence of the clustering of BOSS quasars, we divide the quasar sample into three absolute magnitude subsamples 

covering the same redshift range (2.2 ^ 2 ■. 2.8). The three panels display the real space-correlation function fitted by a one-parameter power law over 

scales of 4 < r p < 25 h ~ 1 Mpc. The red dashed line is the one-parameter fitted power-law with 7 = 2.0 and the blue solid line is the one-parameter fitted 
power-law to the measured correlation function of the main sample plotted here as a reference line. 
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Figure 15. Similar to Fig.[To]but for samples in the redshift range 2.2 $1 z ^ 2.8, subdivided by absolute magnitude. See Table[5]for the fitting results. The 
green curve is the fitted dark matter correlation function to the measured £( s) for each sample and the overlaid blue dashed curve is the fitted curve to the §(s) 
for the main sample, plotted here as a reference for comparison. 


to larger scales than we are willing to trust given our observational 
syst ematics (s < 25 ft -1 Mpc), and it is possible that the lShen et alJ 
d2007l) results on larger scales are somewhat contaminated by sim¬ 
ilar systematics that are obfuscated by the larger statistical errors 
associated with their measurement. With this_caveat in mind, how¬ 
ever, we will take the results of lShen et alJ j2007l) at face value and 
assume that their measured clu stering on l arge s cales is accurate. 
If the high biases measured bv lShen et all 120071) are correct, they 
would im- ply a sharp change in host halo masses at luminosities 
or redshifts slightly beyond the boundaries of the samples analyzed 
here. 


To compute the duty cycles in Table[7] we compare the cumu¬ 
lative luminosity function of quasars over a range of luminosities 
to the c umulative space dens i ty of halos over a range o f host halo 
masses jHaiman & Huill200ll : lMartini & Weinberg|200ll) : 


/duty — 


f r Lm “ 5>(L)dL 

J ^min_ 

f 

JI 


M h ,, 


r dM 


( 10 ) 


We take the quasar luminosity function < D(L) from lRoss et alJ 
l i2013ai) . and report values of the number density based on this 
luminosity function in Table [7] iRoss et all d2013al) derive careful 
corrections to the number density of BOSS quasars by simulating 
quasars and attempting to reselect them with the BOSS targeting 
algorithm. This work on the BOSS selection function implies that 
BOSS only recovers ~ 65% of 2 ~ 2.5 quasars at i ~ 20 falling 
to o nly ~ 35% o f z ~ 2.5 quasars at i ~ 21. (see Table 4 of 
iRoss et aT]|2013al) . This selection function is reflected in the dif¬ 
ferences in our characteristic number densities in Table [6] and the 
number densities derived from d?(L) in Tablc[7] In concert, the val¬ 
ues from these two tables imply that our highest luminosity cluster- 
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Figure 16. The variation of the measured correlation length and measured and expected bias as a function of bolometric luminosity for quasars in three 
absolute magnitude subsamples. Each of the luminosity subsamples has a similar redshift (z = 2.456, z = 2.436 and z = 2.411). The number density of 
quasars in each sample (n) is recorded in the legend of the left-hand panel. The blue solid line is the best-fit power law (in the general form of a(L/Lo) m 
with logio (Lo) = 46.33) to the measured correlation length and bias values in the three luminosity bins. The black dashed line is the best-fit power law with 
m — 0 (i.e. a model for which clustering does not depend on luminosity). In the right-hand panel, the coloured symbols are the best fit biases to the measured 
£(s) for each sample. The slopes of the fitted line to the measured (coloured) biases and correlation lengths are m = —0.027 =b 0.27 and m = 0.56 =b 0.80 
respectively. The connected grey points in the right hand panel show the predicted luminosity dependence of bias if quasar luminosity is monotonically related 
to halo mass and all halos have the same quasar duty cycle / = 0.0073 that best fits the full sample. 


ing sample is only ~ 55% complete and that our lowest luminosity 
clustering sample is only ^ 30% complete, in good agreement with 
Rosj_et all l i2013;l) . 

Note that it is reasonable to integrate our halo masses to 
A7h,max = oo in Eqn. [TO] rather than to some other threshold 
value of Mh,ma.x- This is because, in the formalism of Eqn. [TO] 
we are making no strong assumption that quasar clustering corre¬ 
lates strongly with quasar luminosity. In effect, we are assuming 
that quasars of a specific luminosity could occupy a broad range 
of halo masses. We relax this assumption later in this section by 
inverting our argument in order to predict the biases that would 
arise if halo mass was a tight, monotonic function of luminosity. 
In Tabic [7] we show our results for /duty, derived from Eqn. 1101 
for our main sample and for each of our redshift- and luminosity- 
based subsamples. Given the nearly constant values of Mh^min we 
derive from clustering, one can interpret the /duty values as the 
probability that a halo above Mh ~ 10 12 Ii~ 1 Mq hosts a quasar in 
the corresponding luminosity range, and these probabilities are ap¬ 
proximately proportional to the space densities of each subsample. 

We derive error bars on /duty by matching the ±1<t values of 
ATh.min inferred from the uncertainties. Figure [18] plots our de¬ 
rived values of /duty against redshift for our primary sample and 
for our samples binned by redshift. We do not plot z < 2 values 
because the implicit assumption of a monotonic relation between 
luminosity and halo mass is likely to become inaccurate at red- 
shifts below z k, 2 where the quasar luminosity function begins 
to decline. Our bias measurements imply a duty cycle of approx¬ 
imately 0.7% for our full sample, with a lower value of 0.1-0.2% 
for our highest redshift subsample. Disentangling whether this re¬ 
sult is driven by redshift or luminosity is difficult. By virtue of the 
magnitude-limited nature of BOSS our highest redshift subsample 
also contains moderately more luminous quasars. It is plausible that 


the most luminous quasars should have smaller dut y cyc les. Our re¬ 
sults are reasonably consistent with Fig. 14 of lWhite et alj d2012h . 
which suggests that duty cycle is only weakly dependent on lumi¬ 
nosity at the level of 0.1-1%. However, we caution that the transla¬ 
tion from quasar bias to duty cycle relies on the assumption that the 
sharp threshold in quasar luminosity corresponds to a sharp thresh¬ 
old in halo mass. 

We can again compare our results to those of Ishen et all 
d2007i) . which imply higher du ty _ cycl es of ~ 10% a t 

2 = 3 dWhite, Martin i & Cohnl [20081 : IShankar et alj l201dh . 

White. Martini & Cohnl d2008h quote 6 q ~ 7 or 6 q ~ 8 for the 
Shen et alj d2007l ) results usin g ‘all’ data or just the ‘best' data. 
In our formalism, taking the Ishen et al. d2007h value of 4> = 
5.592 x 10 -7 (fo^ Mpc) -3 at z ~ 3 (their Table 6) and adopting 
our mass function dTinker et alj|2008i ) and cosmology, we derive 
/duty ~ 2% for 6 q ~ 7 and /duty ~ 6% for 6 q ~ 8. The factor of 

3 difference in these duty cycles is driven by the steep relationship 
between bias and mass as halos become increasingly rare in the 
iTinker et alj d20ld ) model. In contrast, we measure /duty ~ 0.2% 
at 2 ~ 3 for our sample. 

The broad reason for this order-of-magnitude difference is that 
ishen et all d2007h measure far stronger clustering than we measure, 
implying that their quasar sample is far more biased than our sam¬ 
ple. Their higher bias implies that their quasar sample occupies 
more massive and therefore far rarer halos than our sample. Al¬ 
though their quasars are brighter, and so less numerous than our 
quasars by a factor of 7, their measured bias implies that their 
quasars are in halos that are a factor of ~ 70 rarer than our ha¬ 
los (for their ~ 7 re sult). In other words, the lower density of 
the brighter Ishen et alj (20071 ) quasar sample compared to our own 
does not compensate for the far lower density of the occupied ha¬ 
los implied by the Ishen et alj d2007i ) measured bias. Taken at face 




















































16 


Eftekharzadeh et al. 


CU 

O 



□ Main sample 

o Redshift subsample 1 
A Redshift subsample 2 
V Redshift subsample 3 
o Font-Ribera+2013 

□ Shen+2007(good fields) 

□ Shen+2007(all fields) 


* ® $ 


!! 

[] 


i 

i 


■. [ i 

ii- 


Filled symbols: M h 
Open symbols: M hmin 


2.25 2.83 _ 3.57 4.50 

z 


Figure 17. The evolution of minimum and characteristic halo mass [defined respectively by bo = b( M > Mh, m in) and 6q = b(Mh)] for the main samp le 
and three r edshift subsamples (deta iled in Table [4} using the halo mass function from Tinker et al] J2010I) . We also measure these quantities for lshen et afl 
I 2007 I) andlFont-Ribera et al.l 120131) based on their reported bias values. See Table[7]for the values and their 1-sigma uncertainties. The ’’good fields” points 
have been offset slightly for visual clarity. 


value, the combination of our results with Ishen et all (120071) im¬ 
plies a sharp drop in the quasar duty cycle near 2 « 3. As we 
have previously noted in this section, other possible reasons for 
this discrepancy might include the difficulty in measuring quasar 
clustering for faint, high-redshift quasars, as well as uncertain- 
tie-in the theoretical halo mass and bias function at high redshift 
i jlinker et alfeoidi . 

If quasar luminosity is a monotonic function of halo mass, 
then more luminous quasars should be more strongly clustered at 
a given redshift. As we have discussed previously in this section, 
this seems unlikely given our measurement that bias is not a strong 
function of luminosity. We can determine the biases that would be 
implied by a monotonic relationship between halo mass and lumi¬ 
nosity by fixing /duty in Eqn. [TO] We can then assume that our 
brightest quasar sample occupies a mass range of Mh , m in < M < 
00 , integrate over this range to determine Mh,m in and then use this 
Mh, m i n as the maximum halo mass Mh,max for our next most lu¬ 
minous sample. In this manner, we can derive Mh ranges for each 
of our luminosity subsamples and then integrate over those ranges 
(from Mh,m in to Mh, max rather than always from Mh,m in to 00 ) 
in Eqn.[9]to derive bias values. 

To begin, we assume that the host halos for each of our lumi¬ 
nosity subsamples at z = 2.4 have the value of /duty = 0.0073 
that we derived for the full sample at this redshift. For our high¬ 
est luminosity bin, we then solve for the value of Mh,m in from 
Eqn. 1 101 setting L m i n = 2.99 x 10 46 ergs -1 (see Table [7} and 
/duty = 0.0073, obtaining Mh,m in = 2.39 x 10 12 /i _1 Mq. For 
the next bin, we use Eqn. GSI With /duty = 0.0073 and L min — 

1.39 x 10 46 ergs _1 , but now we set the upper limits of the in¬ 
tegrals to the lower limits L m i n and Mh, min for the higher lumi¬ 
nosity bin, deriving a halo mass range Mh/10 12 /i _1 Mg = 1.39- 

2.39 for this luminosity bin. We repeat the process for the low¬ 
est luminosity bin, with L m ax and Mh, max equal to the minimum 


values for the middle luminosity bin, deriving a halo mass range 
Mh/lO 12 /i _1 M 0 = 0.99-1.39. 

These calculations of Mh,m m> where we impose a duty cy¬ 
cle and strictly require that more luminous quasars reside in higher 
mass halos, are conceptually different from those reported in Table 
[7] where we derived Mh, m i n for each luminosity sample indepen¬ 
dently based on its measured clustering bias. Here we are deriving 
Mh, m in and Mh, max from the luminosity function, having chosen 
a duty cycle that will lead to the correct average bias across the 
full sample. Having derived Mh, min and Mh, max we then compute 
the predicted bias for each of our subsamples binned by luminosity 
using Eqn. [9] 


We show the three predicted bias values as connected points 
in the right-hand panel of Fig. m Previous quasar clustering 
measurements at lo wer redshifts have shown little or no tren d 
with luminosity (e.g^ ShentdalJ 200Sk iPadmanabhan et al.ll2009h . 


IShankar, Weinberg & Shenld201C ) found that these results required 
a tight relation between luminosity and mass for low values of the 
duty cycle (~ 0.0001) but quite considerable scatter in the mass- 
luminosity relation for higher values of the duty cycle (~ 0.001). 
In essence, the meaning of the duty cycle itself therefore becomes 
complex in the presence of substantial luminosity scatter (see also 
[Shankar, Weinberg & Miralda-Escudel l2013lf ). Further, we should 
be cautious about the exact value of /duty inferred for our full 
sample, which is used to derive the results in Fig. M because it 
is likely that some lower mass halos host some of the quasars in 
our sample. We must therefore be careful not to strictly interpret 
our results in Fig. |T6] as completely excluding a monotonic rela¬ 
tionship between bias and luminosity. With these caveats in mind, 
however, the high value of we have measured for quasar cluster¬ 
ing near z ~ 2.5 implies that we are on a steeply rising portion of 
the b(Mh) relation, and the small statistical errors in Fig.ll6lsug- 
gest a near-constant value for b(L) and a correspondingly constant 
Mh,min(L). Taken at face value, then, our results imply that all ha- 





























Clustering of intermediate redshift quasars 17 


0.010 

S' 


0.001 

2.2 2.4 2.6 2.8 3.0 


Figure 18. The evolution of the duty cycle for the main sample and three 
redshift subsamples (detailed in Table[4} using the halo mass fu nction from 
iTinker et alJ J2010h and the quasar luminosity function from iRoss et al] 
l2013al) . We also measure /duty for iFont-Ribera et al.l J201 3l) . See Table 
[7]for the /duty values and their lcr uncertainties. 
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los above Mh ~ 10 12 are able to host a quasar anywhere 

in our sample luminosity range, with little correlation between lu¬ 
minosity and halo mass. 

A number of more sophisticated approaches to quasar clus¬ 
tering are poss ible, including abundance matching to galaxies 
dConrov & Whitej 20Jj|l. halo oc cupation distribution ( HOD) mod¬ 
eling (e.g. Degraf et alJ 1201 lbl) . light curve models jLidz et alJ 
120061) . models that add scatter to the luminosity-halo mass rela¬ 
tion dshankar, Weinberg & Shenjl20ld) . and detailed evolutionary 
models t hat track black hole growth and model their re lation to ha¬ 
los (e .g. Ishankar et alJl2O10l) . We refer the reader to I White et alJ 
J2012I) for a more in-depth discussion of some of these issues. The 
quasar bias factor, while a single number, provides a critical quan¬ 
titative constraint for any of these models. We have now measured 
this constraint with high precision at several luminosity thresholds 
over the redshift range 2.2 ^ 2 ^ 3.4. 


7 CONCLUSIONS 

We measured the real and redshift-space correlation functions for 
a uniformly selected sample of 55,826 NGC-CORE quasars from 
the final data release of SDSS-IIl/BOSS over the redshift range 

2.2 ^ z 2.8. We also investigated the luminosity dependence of 
quasar clustering by splitting this NGC-CORE sample into three bins 
of absolute magnitude containing approximately equal numbers of 
quasars. We do not detect a significant luminosity dependence to 
clustering strength for BOSS quasars over a factor ~ 10 range in 
luminosity. Our clustering results are summarized in Table[5] 

This work provides precise quasar clustering measurements 
near z ~ 2.5, using samples of tens-of-thousands of quasars. 
We more than double the sam ple size compared to the earlier 
BOSS work of lWhite et alJ d2012h . and extend measurements of the 
real-space projected 2PCF to larger scales. Our best-fit correlation 
lengths and slopes are in good agreement w ith, but more accurate 
and more precise than, those of lWhite et al.l J20121) . 

In order to investigate the redshift dependence of quasar clus¬ 
tering, we extended the upper redshift limit of our BOSS sample 
to z — 3.4 and divided the resulting set of 73,884 quasars into 
three redshift subsamples of roughly equal size. Our fitting results 
for the three redshift subsamples (see Table [5}, suggest that the 


correlation length does not evolve strongly over the redshift range 

2.2 ^ z ^ 3.4. When compared with earlier work (see Fig. EU, 
this result suggests that ro evolves by no more than a factor ~ 2 
from z ~ 3 to z ~ 1. We fit our redshift-space correlation func¬ 
tions to a p ower spectrum for dark matter halos from HALOFIT 
(Smith et al.l 120031 ) modified by the correction for redshift space 
distortions from Kaiser (1987). The bias of our quasars does not 
depend strongly on either redshift or luminosity within our sample 
and is roughly &q ~ 3.5 at z ~ 2.5. 

We adopted the halo bias model of lTinker et al.l J2010I) in order 
to use our measured biases to estimate the average halo masses of 
quasars in our main (NGC-CORE) sample and in our three redshift 
subsamples (see Table [7}. The estimated characteristic dark mat¬ 
ter halo mass that we derived for quasars over the redshift range 

2.2 ^ z < 3.4 covers a relatively wide range of halo masses from 
~ 0.6 to 3 x 10 12 /i _1 Mq because a higher mass is required to 
yield the same bias at a lower redshift. By integrating o ur adopted 
bias model over the halo mass function of lTinker et al Jl (200 8i) we 
determined the minimum halo mass corresponding to our measured 
biases, deriving Mh , m in values from ~ 0.3 to 1.5 x 10 12 /i -1 Mg. 

By comparing the number density of halos above this mass 
threshold from [Tink er et al.| (20081) to the quasar luminosity func¬ 
tion for BOSS from Ross et al. l 2013al) we found that the duty cy¬ 
cles for our main and redshift subsamples are /duty ~ 1%, with 
our highest redshift (and therefore most luminous) subsample ap¬ 
proaching /duty ~ 0.2%. These values of /duty ~ 1% are in agree¬ 
ment with previous works, and are broadly consistent with q uasars 
being short-lived and episodic (as argued by, e.g., ICroom et alJ 
120051) but are also not inconsistent with a mixed population of 
quasars—consisting of a bursty population igniting as their par¬ 
ent halos grow, combined with a l arger fractio n of shorter-lived 
quasars (see, e.g., the discussion in [White et al.[|2012l) . These val¬ 
ues of /duty implicitly assume that halos with mass below Mh , m in 
do not host quasars in our sample, and they then describe the prob¬ 
ability that a randomly selected halo above Mh , m in is active above 
our sample luminosity threshold at any given time. 

There are two somewhat surprising results in our analysis. 
The first is a nearly constant strength of clustering with luminosity. 
We investigated a model in which the quasar luminosity is a 
monotonic function of host halo mass, using the IRoss et all 
J2013b!) luminosity function to infer halo mass ranges from quasar 
space densities assuming the duty cycle (/duty ~ 0.7%) derived 
from fitting the full sample. This model predicts a quasar bias 
that rises from ~ 3.1 for our faintest sample to 6 q ~ 4.5 for 
our brightest sample, clearly inconsistent with our measurements. 
We conclude that there must be substantial scatter between 
halo mass and quasar luminosity in this redshift/luminosity 
range. A broad short-timescale scatter in Eddington ratio has 
been proposed both from observat ional and theoretical stand 


Kauffmann & Heckman 2009; 

Novak, Ostriker & Ciotti 

2011; 

Aird et al. 2012; Bongiorno et al. 2012; Gabor & Bournaud 

20li Hickox et alJ 20l4 

Veale, White & Conrov 2014); 


hawinski et all [20150 , which further argues against a tight, 
monotonic relationship between luminosity and mass. Our conclu¬ 
sion is stronger than that reached from previous studies at lower 
redshift because of our high statistical precision and because the 
high bias of our quasars puts us on the steeply rising portion of the 
b(Mh ) relation. 

The second surprise is the difference between our measured 
host halo masses and duty cycles when compared to the work of 
IShen et al.l (20071) . At z ~ 3 we derive halo masses and duty cy- 
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Table 7. The first three columns are th e cha racteristics of each subsample; the redshift range, the luminosity range and the number density of quasars calculated 
using the luminosit y function of lRoss et alJl2013ah . The 4th and 5th columns are the minimum and the characteristic halo mass based on the halo mass function 
< Tinker et all2008l) and bias model iTinker et al.l20ld) at the average redshift of the sample (see Eqn.[9j. The 6 th column lists the duty cycle for the quasars 
in each sample, which is derived from Mh , m in and <f> (see Eqn. ED. /duty is expressed as a fraction of the Hubble time (9.785 h 1 Gyrs). 


cles that are close to an order-of-magnitude smaller th an those of 
IShenet aD d2007l) . Both of these results are driven by IShen et alj 
(2007) measuring that quasars cluster far more strongly than we 
have found in this work. The high biases found by IShen et al J 
( 2007 ) imply that the quasars they use are in far rarer halos than 
thos e we study, by a factor of ~ 70 x. The quasars studied by 
[Shen et aO 1200 7) are brighter than the BOSS sample we study, and 
so they are also r arer, but only by a factor of ~ 7 x. The fact that the 
IShen et alj d2007h quasar sample has a number density nearly 10 x 
lower than our BOSS sample but an implied halo abundance nearly 
100 x lower is the main discrepancy between this earlier work and 
our own results. A strong dependence of quasar clustering on lu¬ 
minosity or redshift could explain the difference of results, but no 
hints of such dependence are seen within our sample. 


Perhaps the most likely explanation of these differences lies in 
the difficulty of studying clustering of high redshift quasars targeted 
near the limits of an imaging survey. This is particularly true near 
2 ~ 3, where quasar colours resemble those of stars. Our analy¬ 
sis takes advantage of improvements in SDSS photometric calibra¬ 
tion and target selection algorithms as well as much larger numbers 
that afford greater measurement precision. The excellent agree¬ 
ment between our measured quasar bias and the value found by 
iFont-Ribera et alj d2013l) from cross-correlation with the Lyman- 
a forest is strong evidence that our full sample result, at least, is 
not affected by variations in target selection efficiency or redshift 
completeness (which would both affect auto-correlations but not 
cross-correlations). Nonetheless, the difference from earlier results 
at 2 > 3 highlights the importance of studying large quasar sam¬ 
ples selected from deep imaging with well controlled systematics. 


This paper reports measurements of quasar clustering using 
the final sample from the SDSS-III/BOSS survey. The prospects for 
further refining measurements of quasar clustering across a range of 
redshifts are excellent. The SDSS-IV/eBOSS survey has just begun 
(see Myers et al. 2015) and will ultimately result in a uniformly 
selected sample of ~500,000 quasars over a range of redshift from 
2 ~ 0.7 to 2 ~ 3.5. SDSS-IV/eBOSS should more than double the 
number of uniformly selected quasars at moderate redshift (2.2 < 
2 < 3.5) compared to the sample used in this paper, and it will 


provide an opportunity to measure quasar clustering across most of 
cosmic history in a single, consistently selected sample. 
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APPENDIX A: THE “BOSSQSOMASK” SOFTWARE 

In order to make clustering measurements it is necessary to cre¬ 
ate a random catalog that mimics the selection function of the 
data but that is otherwise unclustered (see Sj2}. Here we introduce 
the “BOSSQSOMASK” package—a series of reasonably mature IDL 
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codes that can be used to produce random catalogs that mimic the 
BOSS CORE quasar selection. BOSSQSOMASK has been used in pre¬ 
vious papers fe.g.. Iwhite et alJl2Q12l : fKaragiannis. Shanks & Rossi 
12014 but we make the current version available onlin^l with the 
publication of this paper. Each BOSSQSOMASK routine has a header 
documenting its purpose, but here we highlight the critical steps 
needed to produce coordinates and redshifts for random objects 
according to the angular spectroscopic completeness and redshift 
distribution of BOSS quasars. 

To mimic the selection of BOSS quasars, approximately 
300,000 distinct spheric al caps are created using the MANGLH 11 ! 
software dSwanson et al.l 120081) . To cull areas of the survey that 
greatly affect the quasar density, we use veto masks developed 
to analyze BOSS galaxy clusterinj^l (se e, e.g.. Iwhite et al.ll201 ll : 
I Anderson et alj 120121 : iRoss et all l2012al> . The angular complete¬ 
ness of the survey is tracked in regions defined by a unique set of 
overlapping spectroscopic tiles (called a sector; see iBlanton et alj 
l2003i) . The angular completeness in each sector (called / comp in 
this paper) is defined to be the ratio of the number of quasars for 
which good spectra are obtained to the number of all targets in the 
target list. The BOSSQSOMASK code takes the minimum accept¬ 
able angular completeness as an input (e.g., / com p=0.75 can be sent 
to indicate that the survey should be limited to areas with angular 
completeness greater than 75%, as was used in this paper). 

In the rest of this appendix, we detail the specific requirements 
and procedures for the BOSSQSOMASK software. There are numer¬ 
ous observational effects that have a sec ond-order effect on target 
density in BOSS (e.g.. IRoss et alj|2012al) . Such effects are particu¬ 
larly prevalent in the SGC region of BOSS (see the eBOSS target se¬ 
lection papers, Myers et al., Prakash et al. 2015). BOSSQSOMASK 
does not attempt to correct for these higher-order effects, and so 
users should not trust the software to account for subtle differ¬ 
ences in target density between the random catalog and the quasar 
sample—differences that can cause deviations from zero clustering 
amplitude on very large scales. 

Al File and system requirements 

BOSSQSOMASK requires a list of all objects that would have been 
targeted as a BOSS CORE quasar, i.e., all point sources in the BOSS 
area that have an XDQSO probability above 0.424 to the magnitude 
limit of BOSS. An appropriate filq3 is included in the BOSSQSO¬ 
MASK data directory, and is used to make the psuedo-CORE of all 
quasars that “should have been targeted” (i.e. if angular complete¬ 
ness was 100%). 

The package requires two environment variables to be 
definec£3: $BOSSQSOMASK_DIR and $BOSSLSS_DIR. These 
variables are necessary so that BOSSQSOMASK can find masks as 
well as code functions. $BOSSLSS_DIR must contain all of the 
masking files listed as part of the galaxy clustering package at, e.g.. 


the binaries have to be properly compiled such that, e.g., man¬ 
gle’s ransack procedure in $IDLUTILS_DIR/src is functional. 
Running the sequence of the following codes will produce a 
random catalog (with the user’s desired size and completeness 
level) using the survey masks and list of CORE targets. 


A2 Main steps 

1) Making the XD CORE target list (make_xd_core . pro): 

This code takes the file of all XDQSO CORE targets as 
an input and constructs the pseudo-CORE target list (called 
data/xdcore/xdcore . fits). The code uses veto masks for 
bright stars, center posts, bad u-columns and bad-photometry 
fields (many of these masks reside in the directory pointed to by 
$BOSSLSS_DIR). 

2) Making the XD spAll file (make_xdspall .pro): 

This step requires the latest version of the BOSS spAll file 
that contains the spectroscopic information for BOSS targets (e.g. 
spAll-v5_7_0). The spAll file can be downloaded by running 
pro/get/get_spall.pro, which limits the spAll file to just 
the relevant information (the “mini-spAll” file) and places it in 
the directory data/spall. 

Once the mini-spAll file is in the correct path, 
make_xdspall.pro cross-matches observed objects in the 
spAll file with targets from the xdcore.fits file (produced 
in the first step), and makes a file of matched targets and spectra 

called data/compfiles/xdspall. fits. 

3) Assigning weights to the sectors of the survey mask 
(make_sector_completeness .pro): 

In order to construct the random catalog with the correct an¬ 
gular density and to select only sectors at the user-specified com¬ 
pleteness level requires a weight that tracks the angular complete¬ 
ness for quasars in each sector of BOSS. This code takes the files 
of targets and spectra, obtains the BOSS sector geometry from the 
$BOSSLSS_DIR and calculates the angular completeness in each 
sector. 

Several local files are written that are more manage¬ 
able than the comprehensive target files. Most notably the 
file data/compfiles/xdspallmask . fits, which 
indicates which sector of the mask each target occu¬ 
pies. In addition, smaller MANGLE polygon files are writ¬ 
ten that contain the subset of sectors that are relevant to 
BOSS quasars (data/compfiles/bosspoly. fits and 
data/ compf iles/bosspoly .ply). 


http://www.sdss3.org/dr9/tutorials/Iss_galaxy 

A working copy of the IDLUTILS software suit£3 is needed, and 

10 http://faraday.uwyo.edu/-admyers/bossqsomask 

11 MANGLE is fully described and made freely available at 

http://space.mit.edu/-molly/mangle/ 

12 http://www.sdss3.org/drlO/tutorials/lss_galaxy.p 

13 e.g. the file xdcore.pqsomidzgtO . 4_targets . sweeps . fits 
in the data/xdcore/ directory 

14 A detailed README file guides the user through the steps of setting up 
the system to work with the package. 

15 e.g. http : / / www. sdss3 . org/dr8 / software/idlutils . php 


4) Constructing the angular distribution of objects in the random 
catalog (make_angular_random_catalog.pro): 

3 This code reads the file of which sectors are relevant 
for BOSS quasars (e.g., bosspoly .ply) as well as the file 
that indicates the mask sector and weight for each target 
(xdspall. fits) and creates files of targets and spectroscopi¬ 
cally confirmed targets based on the user-specified completeness 
(denoted completeness ), called xcLcomp completeness .fits 
and xdspall_comp completeness .fits. 

The code also determines the sectors that have a higher weight 
than the user-specified completeness level, and then angularly pop¬ 
ulates those sectors (according to their weight) with N times more 
total random points than BOSS quasar targets using mangle’s 
ransack procedure. N is specified as an input by the user. 
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5) Assigning redshifts to objects in the random catalog 
(make_redshift_random_catalog .pro): 

Based on spectroscopically-confirmed quasars in the sec¬ 
tors derived by make_angular_random_catalog.pro this 
code generates a random redshift from the cumulative distri¬ 
bution of quasar redshifts. The code outputs final data and 
random catalogs in fits format of the size and at the com¬ 
pleteness level specified by the user. As with all files gen¬ 
erated during the construction of the random catalog, the fi¬ 
nal output files are stored in the directory data/compfiles. 
The files are called data-compcompleteness . f its and 
randoms.compcorapleteness. f its where completeness is 
the user-specified completeness. 

A3 Examining the produced random catalog 

Finally, the procedure diagnostics . pro is provided to help 
diagnose whether BOSSQSOMASK is functioning as expected, 
diagnostics .pro can conduct four different tests: 

(i) Angular distribution test: the angular distribution of data and 
random objects should match 

(ii) Redshift distribution test: the redshift distribution of data 
and random objects should match 

(iii) Weights versus number of random points test: the num¬ 
ber of random points in each sector should be proportional to 
“completeness x area” for each sector at the user-specified com¬ 
pleteness level 

(iv) Weights versus number of objects test: the number of 
data points in each sector should roughly be proportional to 
“completeness x area” for each sector at the user-specified com¬ 
pleteness level 
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